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Abstract. We report the measurement of first (n = 1) and higher order(n = 2-5) harmonic 
coefficients {v„) of the azimuthal anisotropy in the distribution of the particles produced in 
U-l-U collisions at ^/sNN ~193 GeV, recorded with the STAR detector at RfflC. The differential 
measurement of (u„) is presented as a function of transverse momentum (pt) and centrality. 
We also present Vn measurement in the ultra-central collisions. These data may provide strong 
constraints on the theoretical models of the initial condition in heavy ion collisions and the 
transport properties of the produced medium. 



1. Introduction 

The study of azimuthal anisotropy, based on Fourier coefficients, is widely recognized as an 
important tool to probe the hot, dense matter created in heavy ion collisions [T]. These 
measurements are sensitive to early stage evolution of the system. In a picture of hydrodynamic 
expansion of the system formed in the collisions, these final state momentum anisotropies are 
expected to rise due to the initial state anisotropy driven pressure gradients and subsequent 
interactions of the constituents. Specially the differential measurement of azimuthal anisotropy 
have been found to be sensitive to initial condition, thermalization, the equation of state and 
transport coefficients of the medium in the heavy ion collisions. Azimuthal anisotropy quantified 
by the Fourier coefficient of the particle distribution can be written as 



where pT is the transverse momentum, y is rapidity and (p is the azimuthal angle of each particle 
and ipR is the true reaction plane angle, defined by beam direction and impact parameter vector 
between two colliding nuclei. The sine terms in the Fourier expansion vanish due to the refiection 
symmetry with respect to the reaction plane. The reaction plane is not measurable directly 
a priori, so the Fourier coefficients are determined with respect to the estimated participant 
event planes [T], 



where are the generalized participant event planes at all orders for each event. The 
estimated event plane from detected particles is corrected for event plane resolution. The 




(1) 



Vn = (cosn(0 - "^n))- 




measured Vn with respect to for reasonable event-plane resolutions are closer to the root- 
mean-square values [2] than the mean values. The first harmonic coefficient vi, called directed 
flow, has two components rapidity odd directed flow generated due to bounce off motion of 
the nuclei and rapidity even directed flow due to dipole asymmetry in the initial geometry. 
The directed flow measurement should suppress the the contribution from global momentum 
conservation especially at peripheral collisions along with other non-flow contributions. The 
rapidity odd vi{y) component at 200 GeV Au+Au collision is found very small at mid-rapidity 
and becomes significant at forward rapidity [3j. The rapidity even component vi{y) recently 
reported at 200 GeV Au+Au collisions [1] form STAR collaboration is found to have weak 
rapidity dependence as predicted by some models [Fj. The transverse momentum dependence 
of rapidity event vi{pt) is of interest in the present study. The second harmonic coefficient 
V2, popularly known as elliptic flow have been extensively studied both experimentally and 
theoretically. Recently higher order harmonics have also gained attention both from theory and 
experimental community. These higher order harmonics can provide valuable information about 
the initial state of the colliding system [6l [71 [8] and also provide an natural explanation to 
the ridge phenomena in heavy ion collisions [9]. Measurement of these harmonics in different 
systems and collisions energy helps to further understand the heavy ion physics in general. 

In uranium - uranium (U+U) collisions, there is the potential to produce more extreme 
conditions of excited matter at higher density and/or greater volume than is possible using 
spherical nuclei like gold or lead at the same incident energy [10]. Uranium has quadrupole 
deformed shape. So U+U collisions may offer an opportunity to explore wider range of initial 
eccentricities. The collisions of special interest are the "ideal tip-tip" orientation in which the 
long axes of both deformed nuclei are aligned with the beam axis at zero impact parameter, 
and the "ideal body-body" orientation in which the long axes are both perpendicular to the 
beam axis and parallel to each other at zero impact parameter. The "ideal tip-tip" and "ideal 
body-body" collision events allow to test the prediction of hydrodynamical model(s) by varying 
the transverse particle density at spatial eccentricity similar to central Au+Au collisions. 

We report here the first Vn measurements in U+U collisions from the STAR experiment at 
the Relativistic Heavy Ion Collider (RHIC) at y^sJvTV = 193 GeV. These data were collected 
during RHIC run-XII in the year 2012. These measurements are presented as a function of 
transverse momentum (px) and centrality. 

2. The STAR experiment and Analysis Details 

Data reported in this proceedings were collected in U+U collisions at ^/snn = 193 GeV in the 
year 2012 with a minimum bias trigger and the ultra central events were taken with a dedicated 
central trigger. Main detector used in the present measurements is the Time Projection Chamber 
(TPC) [llj, the primary tracking device at STAR. TPC has full azimuthal coverage and uniform 
acceptance in ±1.0 units of pseudorapidity. The charged particle momenta are measured by 
reconstructing their trajectories through the TPC. The events with the primary collision vertex 
position along the beam direction (Vz) within 30 cm of the center of the detector are selected 
for this analysis. Event vertex is further required to be in the transverse direction within 2.0 
cm from the center of the beam pipe. Centrality classes in U+U collisions at ^/snn = 193 GeV 
are defined using the number of charged particle tracks reconstructed in the Time Projection 
Chamber(TPC) within pseudorapidity |r/| < 0.5 and passing within 3 cm of interaction vertex. 
The uncorrected, i.e. not corrected for acceptance and reconstruction efficiency, multiplicity 
(A'ch) distribution for events with a reconstructed primary vertex is shown in Fig. [TJ 

This analysis was carried out on tracks that had transverse momenta > 0.15 GeV/c, passed 
within 3 cm of the primary vertex, had at least 15 space points in the main TPC acceptance 
{\r]\ < 1.0) and the ratio of the number of actual space points to the maximum possible number 
of space points for that track's trajectory was greater than 0.52. 
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Figure 1. Uncorrected multiplicity distribution with |?7| < 0.5 in U+U collisions at ^/snn 
193 GeV. Events with Net > 10 are selected for the present analysis. 



For any Fourier harmonic n, the event flow vector (Qn) and the event plane angle (^'n) are 
defined by [1], 

~ (3) 



Qn COs{n^n) = Qnx = ^Wi COs{n<j)i), 

i 

Qn sin(n^'„) = Qny = ^Wi sin(n</)j), 

i 

^1 Qny 



tan 



Qr. 



(4) 



(5) 



where sums extend over all particles i used in the event plane calculation, and and Wi are the 
laboratory azimuthal angle and the weight for the i-th particle, respectively. The event plane 
vector is reconstructed from tracks with transverse momentum [pt) up to 2 GeV/c and within 
\r]\ < 1.0. For first harmonics(n=l), the weight is taken as, 



PT 



<Pt> 



<PT > 

and for higher harmonics (n > 1) the weight is equal to the transverse momentum 

Wi = PT- 



(6) 



(7) 



The < Pt > and < Pj^ > represent event averaged quantities. The choice of this weight 
corrects the effect of the momentum conservation and also it cancels out the conventional 
pseudorapidity asymmetric directed flow O [12] . 

We used scalar product method as well as the event plane method to measure the signal. In 
the scalar product method, the particles with 0.5 < rj < 1.0 were assigned to one subevent and 



particles with —1.0 < t] < —0.5 to the other subevent, separated by r] gap of 1.0 units between 
two subevents and at least 0.5 unit with the particle of interest and the subevent to which it 
is correlated. In order to remove the acceptance effects we applied recentering correction [1] 
to the flow vectors. With unit vector defined as Uni = e™''^, the first harmonic coefficient Vn is 
evaluated as, 

vu„ > 0) = <'^7'; <»'";■> (8) 

Mv < 0) = '^"•^'^"■f (9) 

V \^naHnbl 

In the event plane method, the particles with 0.1 < < 1.0 were assigned to one subevent 
and particles with —1.0 < rj < —0.1 to the other subevent, so that subevents were separated 
by T] gap of 0.2 units. Larger pseudorapidity separation was desired but we are limited by the 
event plane resolution. In this method, we first evaluate the event plane angle from the event 
plane vector (Qn) and the event plane angle was flattened using a shifting correction method to 
correct for the detector acceptance effects. The particle correlation with event plane angle 
is given by, 

( p^ {cOSnjcb-^naiV <0))) . . 

' V{cosn{^na{r]<0)-^nb{r]>0)))' 

, ^ Qx {cos n{(t)-^nb{'n > 0))) -,x 

^ v/(cosn(^„a(7?<0)-^„,(r/>0)))- ^ ^ 

3. Result and discussion 

Results are presented with only statistical errors. In these studies, contribution from short 
range correlation such as Bose-Einstein correlations, coulomb interactions are studied using 
a pseudorapidity gap of at least 0.5 units in pseudorapidity between the event vector and the 
particle of interest using scalar product method. Results from both methods, scalar product with 
larger psuedorapidity gap and event plane method with smaller pseudorapidity gap are consistent 
with each other, which suggests that non-flow contribution from short range correlations are 
small. To reduce the effects from high pT particles in the estimation of the event plane, we used 
Pt weight only up to 2 GeV/c. The non-flow contribution from the jets/minijets is not known 
and might be a signiflcant contributor to the systematic uncertainties especially at peripheral 
collisions. 

Figure [2] (a) presents the f„ integrated in transverse momentum px (0.15 < pT < 2GeV/c) 
and pseudorapidity rj (|r/| < 1.0) as a function of centrality. Except for V2, we observe very weak 
centrality dependence. The initial overlap geometry, which makes the dominant contribution to 
the elliptic flow, changes from central to peripheral collisions, giving rise to the strong centrality 
dependence of V2- For vi, centrality dependence is observed only at more peripheral collisions 
(40-80%) where non flow contribution may be significant. The non fiow contribution from the 
momentum conservation effect [12j may hot have perfectly suppressed in peripheral collisions 
in the measurement of vi. The centrality dependence of the geometry fluctuations are not 
yet understood. The weak centrality dependence of the flow coefficients which are believed to 
originate from geometry fiuctuations indicates that geometry fiuctuations weakly depend on 
centrality. 

Figure [2] (b) presents the elliptic flow V2 as a function of centrality compared with a model 
prediction based on glauber based [TS] prediction. Model curve is for 200 GeV U-|- U collisions 
from reference ^14j. For central collisions, experimental data are lower than the prediction. 
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(a) Vn integrated in pt < 2 GeV/c and \rj\ < 1.0 as a (b) V2 as a function of centrality compared to the 
function of centrality prediction based on MC glauber model [14] 



Figure 2. Vn ss a function of centrality for U+U Collisions 



The Vn of all charged hadrons, for n = 1, 2, 3, 4 and 5 as a function of pT at various 
centralities, as shown in Fig. 3. The data for 0-5% centrality are shown in the upper panels, 
and for intermediate centrality (30-40%) in the lower panels. All Vn measurements show an 
increasing trend as a function oi pT- Except for V2, we observe very weak centrality dependence. 
Also shown are the same measurement from Au+Au collisions at 200 GeV for the comparisons. 
We observe the difference in Vn for n=l and n=2 at 0-5% central between U+U and Au+Au 
collisions. This may hint to initial overlap geometry difference in the central collisions between 
two systems Au+Au and U+U collisions. The difference diminishes in higher harmonics and 
more peripheral collisions. 
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Figure 3. f„ measurement at 0-5% central in upper panels and mid central (30-40%) collisions 
in the lower panels as a function of px for U+U collisions at 193 GeV compared with Vn 
measurement for Au+Au collisions at 200 GeV [15] 



We also report measurement of u„ for ultra central U+U collisions. The centrality selection 
is based on reference multiplicity. We subdivide 0-5% central bin into 10 smaller centrality 
bins upto most central 0-0.5% centrality. The most central 0-0.5% collisions should have have 
dominant contribution from tip-tip collisions. In Fig. U] (a), function of centrality for 

ultra central collisions is shown. Other than second harmonic coefficient, the Vn do not change 
in this centrality range. We observe small change for V2 since it still has some contribution from 
initial overlap geometry. This observation suggests that there is still a small contribution from 
body-body collisions. 
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(a) Vn as a function of centrality at 0-5% central (b) «„ as a function of pr at ultra central(0-0.5%) 
collisions collisions 



Figure 4. Vn as a function of centrality and transverse momentum for ultra central collisions 



In Fig. 131(b), the vi signal along with V2, fs, f4 and a function of transverse momentum up 
to ~ 5 GeV/c is shown for ultra central (0-0.5%) collisions. In these most central events we 
find that higher harmonics are also significant in magnitude compared with second harmonics 
in the intermediate transverse momentum. Contributions from higher harmonics should not 
be overlooked interpreting dihadrom correlation data at intermediate transverse momentum. 
These higher harmonics may offer natural explanation to the novel ridge phenomena in heavy 
ion collisions [Gj. Model comparisons with these new data may help us to better understand the 
medium properties. 

4. Summary 

We report the first measurement of azimuthal anisotropy Vn for n = 1-5 as a function of 
transverse momentum pT and centrality in U-I-U collisions at ^/s^n = 193 GeV, recorded with 
the STAR detector at RHIC. Centrality dependence is week for harmonics other than second 
harmonics. For higher harmonics and mid central collisions, u„(U-|-U) is similar to f„(Au-|-Au), 
the difference appears at central collisions for vi and V2- At intermediate pT ranges 3-5GeV/c, 
Vn are comparable to the f 2 signal in ultra central collisions. Model calculation specially at ultra 
central collisions may be useful to constrain the initial condition and transport coefficient of the 
medium produced in heavy ion collisions. 
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